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Abstract 



This paper is the first in a series of three which attempt to resolve the diffi- 
culties that have plagued the NN — irNN problem for the past ten years. The 
problems may be summarized by saying that the current NN — ttNN equa- 
tions cannot fit the experimental data. Various theoretical inconsistencies in 
the current formulation have been pointed out and this work aims to eliminate 
these inconsistencies and so, we hope, produce agreement with experiment. 
This is to be done by using covariant perturbation theory, in which these in- 
consistencies are not present. The covariant perturbation theory is developed 
starting from a model Lagrangian, in order to fix notation and phases. (The 
approach is actually Lagrangian independent but we use a specific Lagrangian 
for concreteness.) It is shown that both old-fashioned "time-ordered" pertur- 
bation theory and the convolution integral of Kvinikhidze and Blankleider 
may be recovered from the covariant perturbation theory when certain ap- 
proximations are made. The connection of these results with the work of 
Klein, Levy, Macke and Kadyshevsky is discussed. Two forthcoming papers 
will pursue this covariant calculation in the NN — ttNN system using the 
model and perturbation scheme developed in this paper and derive fully co- 
variant NN — ttNN equations without the double counting problems present 
in previous covariant equations. 
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I. INTRODUCTION 



Models in nuclear and particle physics are of two main types: (i) Those based directly on 
Quantum Chromodynamics (QCD), in which quarks are the basic degrees of freedom, e.g. 
Chiral Perturbation Theory or the Nambu-Jona-Lasinio model [|3|-[(|; and (ii) Those 
in which mesons and baryons are the basic degrees of freedom, e.g. the Walecka models 
of the nuclear many-body system, which are collectively known as Quantum Hadrodynam- 
ics (QHD) [[F]|§. In principle these two types of models are connected, since the Lagrangian 
for QHD should be obtained from that of QCD upon integration of the quark-gluon degrees 
of freedom. However, as yet this has not been achieved, and one has to resort to models 
of QCD such as chiral bag models || in order to explain the connection between the quark 
and the hadronic degrees of freedom. Since the hadronic picture can in this way be justi- 
fied from the quark picture one would expect that the hadronic models would be adequate 
approximations to the true physical situation. Whether this is the case or not can only be 
determined by comparison of the predictions of the hadronic models with experiment. A 
disagreement between theory and experiment can then be due either to some approximation 
in the calculation using the model, or to the need for explicit quark-gluon degrees of freedom 
in the model. With the current state of computer technology, light nuclear systems are the 
only nuclear systems in which we may perform the calculation of experimentally measured 
quantities while retaining some control over the approximations used. Consequently, only 
in light nuclear systems may we discover whether the models of the second type, in which 
mesons and baryons are the basic degrees of freedom, can accurately describe the prop- 
erties of nuclei. If they cannot, it will be necessary to resort to models of the first type, 
with explicit quark degrees of freedom, in order to obtain an accurate description of nuclear 
properties. Therefore, the interaction of pions and photons with light nuclear systems is a 
testing ground for models in particle and nuclear physics. 

The NN—irNN system is one example of a light nuclear system in which this program for 
the testing of models based on meson-baryon degrees of freedom may be pursued. During the 
last twenty years disagreements between theory and experiment have led to an improvement 
in the calculations based on meson-baryon degrees of freedom. The physics content of 
these models has gradually increased as they attempt to reproduce the experimental data 
ever more closely. To illustrate how experimental results and the need for consistency in 
the theory have together influenced the development of these NN — ttNN models, we will 
briefly review the history of the NN — ttNN equations. This will also allow us to focus 
on some of the problems which led to the current investigation. A detailed analysis of the 
theoretical methods and a comparison of their predictions with experiment may be found in 
the recent book by Garcilazo and Mizutani | ]10| |. 

The suggestion that the one pion exchange potential could be treated within the frame- 
work of the Faddeev equations was first made by Varma who calculated AW scattering 
using a one pion exchange potential and the Faddeev equations [PJ| . In this model, the 



absorption and production of the pion is the result of the tiN amplitude in the P\\ channel 
having a bound state pole with binding energy equal to the pion mass. This model was 
later extended by Afnan and Thomas to include pion production in AW scattering, pion 
absorption in nd scattering and nd elastic scattering ||12|| . In this way a unified formulation 
of the three reactions: 
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N + N ^ N + N, 

N + N ^ir + d, (1.1) 

7T + d —>■ 7T + d 

was achieved for the first time. The problem with this model was that the nucleons were 
not identical: one of the nucleons (labeled the N') was treated as a nN bound state and so 
could emit a pion, while the other was an elementary nucleon, and so could not emit a pion. 

However, the success of the model in predicting the iViV phase shifts for large angular 
momenta, and the ability of the model to describe ird scattering by not only summing the 
multiple scattering series, but also including the effects of real absorption, motivated several 
groups to reformulate the model in such a way that both nucleons could emit a pion. All 
of these formulations were based on a field theory allowing the possibility of an arbitrary 
number of pion emissions and absorptions. However, the constraint of unitarity allowed all 
the groups to truncate the field theory. Some groups performed this truncation using Fesh- 
bach projection operators |I3 JT^|, while others did it through the classification of diagrams 



technique of Taylor ||15| . All groups truncated the theory by neglecting or approximating 
processes involving four or more particles — the so-called one-pion approximation — and so 
obtained a set of linear integral equations, known as the NN — nNN equations, which are 
three-dimensional in momentum space flq-p5|. The equations were then solved and a de- 
tailed comparison made with the extensive experimental data generated in the 1980s. Early 
attempts at extending the Hilbert space to include the mrNN states, i.e. truncating the 
field theory at the next stage, led to a more complicated set of equations which included con- 
tributions from the diagram in Fig. [1] |26|-f28|j . These equations simplify to the NN — nNN 
equations when the diagram in Fig. [I] is excluded |29|] . 

Although these Faddeev-type NN—nNN equations resolved the problem of non-identical 
nucleons by dressing both fermions, the dressing of the two-nucleon propagator was still not 
complete. The truncation of the Hilbert space and the use of time-ordered perturbation 
theory meant that both nucleons could not be dressed at the same time. I.e., in the one- 
pion approximation, the diagram on the right of Fig. was included as part of the dressing 
for the NN propagator, but the diagram on the left of the figure was not — in spite of the 
fact that the two diagrams are of the same order in the coupling constant, and are just 
different relative time-orders of the same process. This incompleteness of the dressing, first 



pointed out by Sauer et al. pCfl , gives rise to a weakening of the ttNN coupling constant as 
the number of nucleons in the system increases. In fact, Sauer et al. showed that in models 
which, like the NN — nNN equations, use an incomplete dressing of the multi-nucleon 
propagator, the omission of some of the dressing means that the wNN coupling constant 
goes to zero in the nuclear matter limit. Furthermore, the omission of this dressing also 
leads to a severe underestimation of the pp <->• nd cross-section in theoretical calculations— 
especially in theories which treat the itN amplitude in the Pn channel properly, i.e. as 
the sum of a pole and non-pole term. The only way this underestimation can be remedied 
is to dress both nucleons fully. Recently, Kvinikhidze and Blankleider have addressed the 
question of how to do this and have elegantly demonstrated how the complete dressing can 



be achieved with the use of a convolution integral representation for the NN propagator |3l 

A second associated problem has plagued all calculations based on the iViV — nNN 
equations which use a Pn ttN amplitude that is the sum of a pole and non-pole contribution. 



3 



None of these calculations has been able to accurately predict the tensor polarization, T20, 
for 7r — d scattering. In 1988 Jennings suggested that this discrepancy between theory and 
experiment is a result of the iViV — -rrNN equations not including the diagram in Fig. [I], 
which has since become known as the Jennings mechanism |3^j33[] . Indeed, the mechanism 



depicted in Fig. [I] is merely a different time-order of a diagram which is included in the 
NN — TtNN equations, shown in Fig. |3|. Jennings pointed out that the excluded diagram 
may well cancel part of the effect of the one included in the NN — ttNN equations, thus 
producing the correct result for the tensor polarization in ir — d scattering. 

Since the Jennings mechanism and the extra dressing of the NN propagator are merely 
different time-orders of diagrams already included in the NN — ttNN equations, logically 
they too should be included in any calculation of the NN — nNN system. Otherwise, the 
processes which are occurring on one nucleon restrict the processes which may occur on 
the other nucleon. Further, the failure of the NN — nNN equations to correctly describe 
the experimental data suggests not only that these two omitted mechanisms should now be 
included, but that they must now be included. 

The omission of these two diagrams from the present theory is a result of using a trun- 
cation based on time-ordered perturbation theory and unitarity. One might think that the 
problem could be resolved by truncating the Hilbert Space at two-pion states, since then the 
two diagrams discussed above would both be included. As mentioned above, an example of 



a theory which takes this approach is the model of Stingl and Stelbovics [26428fl. However, 



the equations obtained from this model are computationally difficult to solve. The model 
also suffers from similar problems to those of the current NN — wNN equations, but at the 
three-pion level. One example of this is that diagrams such as those on the left of Figure |] 
are included in the model, whereas diagrams such as those on the right are not included, 
because they involve a three-pion intermediate state. Again, this occurs even though the 
two diagrams are just different time-orders of the same process. 

Clearly then the problem rests not so much with the choice of truncation point as with 
a failure to sum over all time-orders when doing the truncation. Therefore, if we wish to 
include all the relevant physics in our model we must use a perturbation scheme in which 
all time-orders are included automatically. This paper presents such a perturbation scheme 
as follows. 

In Section [TI] the field operators we will be using are introduced and the Lagrangian for the 
theory stated. In Section |T| we explain how to calculate the single-nucleon and single-pion 
Green's functions with this Lagrangian. At this point we introduce an approximation com- 
mon in intermediate-energy nuclear physics — we ignore anti-nucleonic degrees of freedom. 
Once we make this approximation, the theory simplifies considerably. We will show that, 
apart from this approximation, the perturbation scheme developed is a fully covariant theory 
of nucleons and pions as given in Bjorken and Drell ||34fl . In Section fV|we turn our attention 
to the two-nucleon Green's function. It is explained how to evaluate Gnn(p'i, P21 Pij P2) ■> and 
the Fourier Transform of the equal-time Green's function Gnn^' , x[, x^t, xi, X2) is calcu- 
lated. We find that: 

dzdz' 

where Gnn(E' — z',p[, z',P2, E — z.px-.z-.P'i) is calculated using the Feynman Rules for 
Gnn{Pi,P2'>PiiP2)- In Section |V] it is shown that similar results hold for the m —>■ m' 



/dzdz 1 
—— 1 G NN {E' - z',p[, z',])^ E - z,p l7 z,p 2 ), (1.2) 
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particle Green's function, a result of particular relevance to the ttNN system. This use of a 
convolution integral to calculate multi-particle Green's functions bears a strong resemblance 
to the work of Kvinikhidze and Blankleider (KB) |31],[35] and in Section |VT| we explain how 
to derive KB's result from our field theory. Since KB's work uses a convolution integral 
in order to sum a set of time-ordered (or old-fashioned) perturbation theory diagrams it is 
necessary to first derive time-ordered perturbation theory from our covariant perturbation 
scheme. KB's result then becomes a consequence of Equation ( |1.2| ) and the relationship of 
time-ordered perturbation theory to the covariant perturbation scheme. Once this result is 
established, the approximations used by KB in order to derive equations for the NN — nNN 
system |35[] are critically examined, and are found to have considerable inconsistencies. In 
Section [VI I| we consider how to calculate amplitudes from the Green's function Gnn{E' — 
z',p[, z',P2, E — z,pi, z,f>2) using Lehmann-Symanzik-Zimmermann (LSZ) reduction ||36|| . 
We find that the amplitude generated by a sum of all relative time-orders of time-ordered 
perturbation theory diagrams may be calculated by evaluating the amplitude given by the 
Feynman rules in Sections p| and |IH|. However, as mentioned above, the perturbation scheme 
that generates these Feynman rules is merely an equal-time, no-anti-nucleon approximation 
to the full covariant perturbation theory. Therefore, if we are to include the missing diagrams 
in our theory we need to use some sort of covariant perturbation theory in order to derive 
new four-dimensional NN — nNN equations, even if at certain stages in the calculation we 
make approximations. 

The results described in this article are similar to those obtained by Klein, Levy, Macke 
and Kadyshevsky (KLMK) |37|-|48[1. The essential difference between KLMK's work and 
ours is that the KLMK method is driven by the desire to obtain a three-dimensional integral 
equation, whereas one of the main theses presented in this article is that four- dimensionality 
is an essential feature of any theory which is not going to treat different time-orders of the 
same physical process differently. We conclude our arguments for this thesis in Section |VIII 
by briefly describing the work of KLMK in the context of the calculations performed in the 
previous sections. In particular, we highlight the pieces of physics which are missing from 
their work, due to their use of a three-dimensional equation rather than a four- dimensional 
one. 

Haberzettl has also recently published some work on the the use of a convolution inte- 
gral in iV-body scattering, and has derived four-dimensional integral equations for certain 



scattering problems [49|-[51|]. While his work appears similar to that discussed here it is 
approached from a rather different perspective, being based on cluster dynamics rather than 
field theory. 

In summary then, the inability of the current NN — nNN equations to correctly describe 
the experimental data forces the introduction of new physics into the theory of the iViV — 
nNN system: namely, the consideration of different time-orders of diagrams already included 
in the theory. The only natural way to include this missing physics is to use the full 
covariant perturbation theory developed in this paper. The use of such a perturbation 
theory will, of course, result in four- dimensional integral equations. Any attempt to reduce 
the dimensionality of these equations results in a loss of physics, as in the work of KB 
and KLMK. Therefore, in later papers we pursue the use of the full covariant perturbation 
theory and so derive covariant equations for the NN — nNN system. In this paper we 
wish only to concentrate on the perturbation scheme we shall be using, and to establish the 
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approximations in which it reduces to the convolution integral approach, the time-ordered 
perturbation theory and the work of KB and KLMK. 



II. THE LAGRANGIAN 

In this section we define the field theory we shall be using in order to fix the notation 
and normalization used throughout the paper. We define nucleon and meson creation and 
annihilation operators, state the commutation relations they obey and write down the La- 
grangian and hence the Hamiltonian. The Lagrangian we use is similar to that used in 
Quantum Hadrodynamics calculations |7]f|. Note that although we use a specific Hamilto- 
nian the approach pursued in this paper could equally well be applied to any meson-baryon 
Hamiltonian with vertices which allow the baryons to emit mesons. 

Consider second-quantized nucleon field operators ip(x) and ip(x). We can, of course, 
Fourier decompose these Schrodinger representation operators, as done by, for example, 
Itzykson and Zuber and when we do we obtain: 



m = J (Ifs ^ T,[b a (q)u a (q)e^ + dl{q)v a {q)e^% (2.1) 

= / (2^E^?^ ( ^ ( ^ )e "^ + d^My**], (2.2) 

where m is the nucleon mass; E^j{q) is the energy of a nucleon of momentum q and a is a 
collective label for spin and isospin. Here b^(q) and b(q) are the creation and annihilation 
operators for a nucleon of momentum q while eft and d are the corresponding anti-nucleon 
operators and u(q) and v (q) are the positive and negative energy spinors corresponding to 



momentum q, normalized as in Itzykson and Zuber [52], i.e. such that: 



u a (q)u a r(q) = 5 aa >. (2.3) 

Note that E^(q) can be chosen to give non-relativistic, semi-relativistic or fully relativistic 
kinematics. E.g., for full relativistic kinematics choose: 



E N (q) = \Jm? + g 2 . (2.4) 

Since the nucleons are fermions their operators obey the standard anti-commutation rela- 
tions: 

{b a (q), bl(q')} = (2nf^±S^(q-q>)5 aa , } (2.5) 

m 

{d a (q),dl(q')} = (27r) 3 ^±5( 3 \q-q')5 aa ,, (2.6) 

m 

with all other anti-commutators vanishing. 

Similarly, we define meson field operators <pi{x), with i an isospin label, which we expand 

via: 

Ux) = j - -ff -=r [aifae*** + a\{k)e-^)) , (2.7) 
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where oJ n (k) is the energy of a pion of momentum k, which in the relativistic kinematics is: 



u v (js) = si ml + fc2 ; (2.8) 
and a,i(k) (a](k)) destroys (creates) a pion with isospin label % and momentum k. Again, 



this is in accordance with the definitions used by Itzykson and Zuber ||52|| . The physical 
pion operators are then defined by: 

0± = T^(01±^2), (2.9) 

fa = fa- (2.10) 

Note that while the operators fa are self-adjoint the operators cf>± are not self-adjoint. Since 
pions are bosons the operators a* and a\ obey the commutation relations: 

[ ai (k),4(k')) = 2u 7T (k)(27r) 3 S^\k - k')5 w , (2.11) 

(2.12) 

with all other commutators being zero. 

These commutation relations are written for the Schrodinger operators. Naturally, they 
remain the same if we change to some other representation. The two other representations 
used in this paper are the Heisenberg representation, which we shall denote by a tilde above 
the operator in question, e.g. ip, fa and the interaction representation, which we denote by 
a superscript /, e.g. if) 1 , fa . Schrodinger representation operators will remain unadorned. 
These three different representations are connected by the standard transformations involv- 
ing the Hamiltonian of the system. 

We define this Hamiltonian via the Lagrangian density, C. We choose: 

£ = C-V + C-4> + C-int, (2-13) 

where: 

C D (x) = #r)(*7^ - m)ifj(x), (2.14) 

= \{dj{x) ■ d^fax) - m 2 J{x) ■ fax)), (2.15) 



£int(x) = ~ig j d 4 x N d 4 x' N (fx^x'^^f^XM) ■ fax n )T(x - x' N , x - x N , x - x w ). (2.16) 

Here we have chosen pseudo-scalar coupling for the ir — jV interaction Lagrangian density 
Cj(x), and r( the form factor for this interaction. We introduce 

this form factor in order to model the finite size of all the particles involved. It will also 
allow us to remove the divergences from the theory. (See Figure |5] for a depiction of this 
vertex. ) 

This Lagrangian is that used in the Walecka model (or Quantum Hadrodynamics-1) f?||| 
except for two differences: 
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1. There are no counter-terms to do the renormalization with. Instead, we include a 
form factor in the pseudo-scalar interaction Lagrangian, thus allowing us to introduce 
a cut-off and so remove the infinities from the theory. 

2. In the QHD Lagrangian there are fields corresponding to one family of vector mesons 
and one neutral scalar meson. The above Lagrangian has only one family of pseudo- 
scalar mesons, the pions. Note, however, that this Lagrangian could easily be extended 
to include vector mesons. 

The use to which this Lagrangian is put will also be different. The QHD Lagrangian is used 
for doing nuclear structure calculations and consequently is always solved as a bound-state 
problem, often in the mean-field approximation. Conversely, we will be using our Lagrangian 
to do scattering theory — in particular, the scattering theory of the NN — nNN system. 

From this Lagrangian, we may define the Hamiltonian density, Tl, by a Legendre trans- 
formation^ 

dd ip ddoip 5d (p 

and hence obtain the Hamiltonian: 

H = J d 3 xH(x) — H K + H int , (2.18) 



where: 



H. 



K 



J d 3 x^{x)T N {x)4j(x) + J d 3 x(f{x)T n (x) ■ ${x), (2.19) 



lint = +ig J d 3 X ( J d 4 XNd' l x' N d' i X 7T ?p(x' N )j 5 T1p(xN) ■ 0(^7r)r(x — x' N , X — Xn, X — Xtt)). 

(2.20) 

Note that we may define an interaction Hamiltonian density 

««(*) = igj*x„ d%dW*b»W**) ' **r)r<* - *.). (2.21) 

The connection between the position space operators of different representations is now: 

$(x) = e iHx °i/j(x)e- iHx ° (2.22) 
tf{x) = e iHKX \{x)e- iHKX \ (2.23) 

From these we may deduce the results for ip. Furthermore, it is clear that the same equations 
will hold if we replace ip by any component of 4>. 



1 The canonical momenta, it, for the field eft, may not be given by the usual relation, tt = due 
to the presence of the dissipative force in the interaction Lagrangian. If this is the case it means 
that the relation between the Lagrangian and the Hamiltonian is more complicated than that given 
m Equation ( p?L7| ). We ignore this difficulty for the present, since we are only really interested in 
the overall final form of the Hamiltonian, which we may assume is given by EqB.(El8D, (|2l9|) and 
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III. THE SINGLE-NUCLEON AND SINGLE-PION GREEN'S FUNCTIONS 



In this section we first give a brief explanation of how to obtain a diagrammatic expansion 
for the single-particle Green's functions, and how to find the Feynman Rules in coordinate 
and momentum space for the calculation of these Green's functions as a perturbation series 
of Feynman diagrams. This is done primarily to fix normalizations and phases. We then look 
at the form taken by the free Green's functions in the low-energy limit, when anti-nucleonic 
contributions to the Green's function can be ignored. 

We define the one-particle Green's functions for nucleons and pions by: 

G N (x',x) = (0\T$(x')l£(x))\0), (3.1) 
G vji tf,x) = (Q\T(fa(x%(x))\Q), (3.2) 

respectively, where T denotes the usual time-ordering operator. 

Consider either of the one-particle Green's functions. By following the standard pro- 
cedure as outlined in e.g. Itzykson and Zuber |^2| , perturbation series for GV(a/,x) and 
G-xjiix', x) may be derived. We may then use Wick's theorem to rewrite each term in these 
series as a sum of products of the free single-particle Green's functions: 

Gf(x',x) = (O\T^ I (x''0(x))\O), (3.3) 
G^(x',x) = (0|r(#(xO#(x))|0). (3.4) 

Consequently we are led to the standard interpretation of the perturbation series as a series 
of Feynman diagrams. Each Feynman diagram corresponds to an analytic expression given 
by the Feynman Rules for the theory. The Feynman Rules in coordinate space for the 
theory of nucleons and pions obtained from the Hamiltonian of the previous Section are the 
standard ones given in Bjorken and Drell, pp. 224-5 [Q, with two changes: (i) A form factor 
T for pion absorption and emission must be added, and (ii) Rule 4 must be replaced by the 
following rule: 

"4. Assign coordinates x and x' to the external points. If dealing with G n ji(x',x) also 
assign isospin indices j and i to the external points. Then, if a nucleon line joins x (or x') 
to an internal point y, associate with it a factor G^\y,x) (G$(x', y)). If a pion line joins 
x (x'), with isospin label to a vertex y, with isospin index k, then include a factor 
Gj${y,x) (G^(x',y)) for that line." 

Note that the change comes about because we are dealing with Green's functions whereas 
Bjorken and Drell dealt with amplitudes. 

Now we may define the Green's functions in momentum space, via: 

G H (p',p) = J d*xd 4 x'e ip ' x 'G H {x',x)e- ipx , (3.5) 
where H is either a nucleon or a pion. By translation invariance: 

Gh{x', x) = Gh(x' — x) (3.6) 

and so it can be shown that: 

G H (p',p) = (2n)W\p'-p)G H (p), (3.7) 
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where: 



G H (p) = J d A x R e vpx *G H {x R ). (3.8) 
where Gh(xr) = Gh{x' — x) is the coordinate space Green's function, the rules for which 



are given in Bjorken and Drell [p4 



Consequently we may find the momentum space Green's function by first using Wick's 
theorem to get an expansion for Gh{x' — x) in terms of the free nucleon and free pion 
Green's functions, and then inserting the Fourier representations for the free Green's func- 
tions Gx \y' — y) and Gjf^{y' — y) and the form factor T(y f , y, y n ). For the free Green's 
functions we have: 

G^{y'-y) = j^e-^'-^\^ (3.9) 

G.fW -V) = J ^e-^G^ik). (3.10) 

For the form factor we define: 

T(q',q,k) = J AdVA^'e^e^rfe'.M,), (3.11) 

and the inverse of this Fourier Transformation may be used to find the Fourier Representation 
oiT(y',y,y^). 

The formulae thus obtained are substituted into the expression for a particular Feynman 
diagram, the coordinate space integrations are performed and finally the Fourier transform 
of the result is taken using Eq. ( |3.8|) . (For more detail on this procedure see Itzykson and 
Zuber pp. 267-8 |52|].) When this is done the following Feynman Rules for the construction 
of Gh(p) are obtained. 



Feynman Rules for the single-particle Green's function in momentum space 

1. Draw all topologically distinct, connected diagrams with one incoming leg and one 
outgoing leg. 

In each diagram: 

2. Assign to the external legs four- momentum p, and to the internal nucleon and pion 
legs four- momenta k\, k^-, ■ ■ ■ , ki, where / is the total number of internal lines. To 
each 7riViV vertex assign an isospin label /i, / 2 , . . . , /„. If the external legs are pion legs 
assign to the incoming (outgoing) pion leg's endpoint an isospin label i (j). 

3. To the nucleon line with four- momentum q, where q equals some kj, j = 1, . . . , J, 
assign a factor G^\q). 

4. To the pion line with four-momentum k, connecting vertices with isospin labels l e , l a , 
assign a factor G n \j e (k). 
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5. To each pion emission vertex assign a factor: 

H2)(^5T0(27r) 4 r(g', q, -k)5^(q' + k - q); 



to each pion absorption vertex assign a factor: 

(-^^)(z 75 r i )(27r) 4 r(g', q, k)8^(q f -q-k); 

to each NN annihilation vertex assign a factor: 

(-i9)(HBriWr(-q ! , q, -k)5^(k - q - q'); 

and to each NN creation vertex assign a factor: 

(^)(2 75 r0(27r) 4 r(g', -q, k)5^(q' + q-k). 

In each case q (q') is the nucleon four-momentum before (after) the vertex (where 
applicable), k is the pion four-momentum and I is the isospin index assigned to the 
vertex. 

6. Include a factor of (-1) for each closed fermion loop. 

7. Integrate over all internal momenta: 

d 4 h d 4 k 2 d 4 k! 



(2tiY (2tt) 4 (2tt) 4 
and sum over all repeated isospin indices. 

It is these rules we shall always work with in the calculations. But, in order to calculate 
with these rules we need to know what the free single-particle propagators G$(q) and 
C?J%) are. It follows from Eqs.Q, Q, Q and (gj) for the free single-nucleon 
Green's function and Eqs.(|2.7|), ( |3.2| ), ( |3.8| ) for the free single-pion Green's function that: 



Gf (,) = — , (3.12) 

7T 

i.e. the free single-particle Green's functions are the full covariant propagators. (See Ap- 
pendix ^ for the details of this calculation.) 



The no-anti-nucleon approximation 

However, the use of these full covariant propagators results, in most cases, in a very 
difficult calculation. In order to simplify the calculation somewhat, a restriction is placed 
upon Green's functions in this work. We make the approximation that the contribution 
made by anti-nucleons to the Green's functions is negligible, an approximation which should 
be true in the low energy limit. Making this approximation means that we force: 
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v(q) -> 0;v(q) -> 0. (3.14) 

We expect this to be true at the energies we are interested in. It immediately implies that: 

GP(x',x) = (O^ 1 (x'0 (x)\0}6(x' - x ), (3.15) 

as one would expect in the absence of anti-nucleons. It is also found that, as a direct 
consequence of the derivation of the single-nucleon Green's function in momentum space, 
we have: 

G " (,) - E N (?) q }-E N (f) ' {3M) 

in this approximation. (See Appendix [A] for the justification of these two statements.) 

This no-anti-nucleon approximation not only affects the free single-nucleon propaga- 
tor, it imposes a definite time-ordering on any process contributing to the single-nucleon 
Green's function. This occurs because the theta-functions in the free single-nucleon prop- 
agators define the order of any two pion emission and absorption times y$ and y' Q . I.e., 
the absence of anti-nucleons determines which of the two times in any free-pion Green's 
function, G n jf(y',y), comes first and which second. So, we have to replace Gjf^{y',y) by 

Gwfi (y' ,y)@(y'o — yo), where y' is the later of the two times. In momentum space this has 
the effect of replacing Eq. ( |3.13| ) by: 

G$(k) = 1 - - %bl ° - ■ (3.17) 
2uj n (k) k£ - u) n (k) 

(Again, a justification of this fact is to be found in Appendix |A].) 

Furthermore, in the approximation in which anti-nucleons cannot be created, it follows 
from the conservation of nucleon number that no processes other than free pion propagation 
can contribute to the single-pion Green's function. Therefore: 

G vji (k) = G$(k). (3.18) 

By contrast, if we attempt to express G^ in terms of G$ we find that G^ obeys a 
Schwinger-Dyson equation. 

So, the no-anti-nucleon approximation reduces the momentum space propagators from 
full covariant propagators to those given in Eqs.( |3.16| ) and ( |3.17| ) and prevents any process 



other than free propagation from contributing to the single-pion Green's function. 



IV. THE TWO-NUCLEON GREEN'S FUNCTION 

In the previous section we explained how to obtain a diagrammatic expansion and Feyn- 
man Rules in both coordinate and momentum space for the single-nucleon and single-pion 
Green's functions. Since our aim is to develop a theory of the iViV — uNN system we should 
next consider the iViV Green's function and its relation to Gn and G^. 

The two-nucleon Green's function is defined by: 
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G NN (x' 1 ,x' 2 ;x 1 ,x 2 ) = (0\T(^ N (x[)^ N (x' 2 )'ijj N (x 1 )'il) N (x 2 ))\0). 



(4.1) 



The procedure for obtaining the Feynman rules for this Green's function is exactly the 
same as that used for the single-particle Green's functions in the previous section. Upon 
implementing that procedure we find that the rules obtained in coordinate space for the two- 
nucleon Green's function are the same as those for the one-nucleon Green' function, with 
two modifications: (i) In the two-nucleon case diagrams must obviously have two incoming 
and two outgoing external nucleon lines, and (ii) Disconnected diagrams are now permitted, 
provided they contain no sub-diagrams which are not ultimately joined to the external lines. 
(For examples of the type of diagrams which are and are not allowed see Figure ||) 

Now we attempt to obtain the Feynman Rules in momentum space for the two-nucleon 
Green's function, but with a restriction that the initial times on both particles are equal, 
as are the final times. It turns out that this restriction of equal times generates a form for 
the Green's function bearing a remarkable similarity to the convolution integral used in the 
work of Kvinikhidze and Blankleider (KB) [|3l],[3f|. The argument we give is valid regardless 
of whether or not we are using the no-anti-nucleon approximation. 

Consider the two-nucleon Green's function Gnn(xi, x 2 \ x[, x' 2 ). Suppose that we restrict 
the times the particles begin and end their propagation, so that the two nucleons begin their 
propagation at equal times, i.e.: 



xl = xl = t, (4.2) 



and end it at equal times, i.e.: 



t'. (4.3) 



Then we denote the Green's function by: 

G NN {t',x[,x' 2 \t,x u x 2 ). (4.4) 

Now consider the Fourier Transform of this Green's function. Because there is only 
one initial and one final time, when the Fourier Transform with respect to energy is taken, 
instead of transforming with respect to the individual particle variables p\, p 2 , Pi' and p 2 \ 
which are conjugate to x\, x 2 , x^ and x 2 , we must transform with respect to variables E 
and E' which are conjugate to t and t'. Therefore, 

G NN (E',p[,p 2 ; E,p u p 2 ) 

= Jdt dt'd 3 x 1 d 3 x 2 d 3 x[d 3 x 2 e i ( E ' t '~^<-^< ) G NN (t', x[, x' 2 \ t, x 1; ^e^-^ 1 ^^ (4.5) 
= J d 4 xi d*x 2 d 4 x[ d 4 x' 2 e^ E '^~^<-^b(x% - x°')GWK, x' 2 ; x 1 ,x 2 )5(x° 1 - x° 2 ) 

xe -i(Ex°-pi-X!-p2-X2) (4 6) 

where we have recalled that x\ = t and x®' = t' and then used the two-time Green's function 
Gnn(x[, x 2 ; xi, x 2 ) in place of the equal-time Green's function. Inserting: 

6(a-b) = J^ lz{a ~ b \ (4-7) 
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and then performing the Fourier transforms over Xi,x 2 ,x' 1 and x 2 gives: 

/dz dz 1 
j—^G NN (E' - z',p[, z',p 2 ; E - z, Pl ,z,p 2 ), (4.8) 

where Gnn(Pi ,Pi,P 2 tP^Px-iPIiP^iP^) 1S the Fourier transform of Gnn(x'i,x' 2 ;xi,x 2 ). At 
this stage it becomes clear that E and E' are the initial and final total energy of the system, 
respectively: 

E = p\+pl (4.9) 
E' = p°; + pt (4.10) 

It can also be seen that the integral over z and z' corresponds to an integral over all possible 
initial and final "relative" energies of the two-nucleon system.^ 

Now the translational invariance of Gnn(x[, x' 2 ; x±, x 2 ) implies momentum conservation 
in G NN (p' 1 ,p' 2 ;p 1 ,p 2 ), i.e. 

G N n(p' 1 ,P 2 ;Pi,P2) = (27T) 4 5 {4) (p' 1 +p 2 ~p 1 - p 2 )G NN (p' 1 ,p 2 ;p 1 ,p 2 ). (4.11) 
This suggests that: 

G NN (E' - z',p[, z',p 2 ; E - z,p!,z,p 2 ) (4.12) 

= (2n) 4 5(E' - E)5^{p' 1 + p' 2 - Pl - p 2 )G NN {E - z' ,p[,z' ^ E - z, Pl ,z,p 2 ). (4.13) 

Further, translational invariance of G^jsf(t', x[, x 2 , t, x\, x 2 ) implies: 

Gnn(E' ,p[,p 2 ; E,p,,p 2 ) = (2tx)%E' - E)5®(p[ + p' 2 - Pl - p 2 )G NN (E,p[, p' 2 , Pl ,p 2 ). 

(4.14) 

These last two results suggest that Eq.( [4.8[ ) may be rewritten as: 

/dz dz' 
-j—-G NN {E - z',p[, z',p 2 ; E - z,px, z,p 2 ). (4.15) 

Therefore the problem of calculating G nn(E , p[, p 2 , p±, p 2 ) is reduced to that of calculating 
Gnn(p'i,P 2 ]Pi,P2)- But Gnn(p'i,P 2 ',Pi,P2) is to be calculated using the Feynman Rules in 
momentum space for the two-nucleon Green's function. To obtain these Feynman rules from 
those in coordinate space we merely proceed as we did above for the one-nucleon Green's 
function. This leads to Feynman Rules for the calculation of Gnn(E — z' ,p[, z' ,p 2 , E — 
z,pi,z,p 2 ) which are exactly the same as those for the single-particle Green's function, 
except that Rules 1 and 2 become: 



2 Note that the definition of "relative" energy used here differs from the standard k° = p® — p 2 . 
Instead, when the term relative energy is used in this paper it refers to the energy z = p 2 = E — p®. 
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1. Draw all topologically distinct diagrams with two incoming and two outgoing legs, ex- 
cluding diagrams in which any sub-diagram is not ultimately connected to an external 
line. 

For each diagram: 

2. Assign to the incoming (outgoing) external nucleon legs energies E — z and z (E — z' 
and z') and momenta p\ and f>2 (p[ and p£). Assign to the internal nucleon and pion 
legs momenta hi, . . . ,kj, where I is the number of internal lines. Note that once the 
momenta and energy of the particles in the initial state are chosen there are two ways 
of choosing the momenta and energy of the particles in the final state. Both possible 
assignments must be included as separate diagrams. The two resulting diagrams differ 
by a factor of (-1), since the nucleons are fermions. 

This is the only change which needs to be made to the full covariant theory in which 
anti-nucleons are included. However, if we are using the no-anti-nucleon approximation 
another change is needed in order to maintain consistency in our work, as follows. Recall 
that, in the calculation of the single-nucleon Green's function we had to replace Gjj} \y' — y) 

by Gnjiiv' ~ y)6(y'o ~ yo), because the no-anti-nucleon approximation imposed a definite 
time-order on pion propagation. However, in the two-nucleon Green's function there is noth- 
ing which specifies the relative time-order of processes taking place on different nucleons. 
Therefore, pions emitted on one nucleon and absorbed on the other do not have the time- 
order of their emission and absorption specified, whereas, as described above, pions emitted 
and absorbed on the same nucleon do have that time-order specified. So, in order to treat 
transmitted pions and pions that are reabsorbed by the same nucleon which emits them 
equivalently, we split the transmitted pion's Green's function G 7r jf(y 2 , y%) into two possi- 
bilities, according to whether y® < y® or y\ < y\, and insert appropriate theta functions. 
Therefore, for transmitted pions one diagram in which we have Gjf)(y', y)is replaced by two 
diagrams, in one of which we have Gjf^{y', y)9(y' — yo), and in the other of which we have 

G-wji (y' \y)0(yo — y'o)- This rule ensures that the relative time-order of a particular pion's 
emission and absorption is always determined. It is important to note that, other than this 
restriction, the relative time-orders of processes on different nucleons remains undetermined. 

If we use the no-anti-nucleon approximation and make this change in order to maintain 
consistency then we still have the same rules for the two-nucleon Green's function as in the 
full covariant theory, but with the simplified propagators given in Eqs.(|3.16|) and (|3.17|). 



V. THE MULTI-PARTICLE GREEN'S FUNCTION 

The argument for the m —>■ m' particle Green's function is constructed exactly as are the 
above arguments for the one and two-nucleon Green's functions. The only change necessary 
to the Feynman rules for the two-nucleon Green's function in coordinate space is that we 
must now draw all possible topologically distinct diagrams with m incoming and m' outgoing 
legs. 

Now we examine the equal-time Green's function, in order to derive a result analogous 



to Equation Q4.8Q . The argument used is similar to that by which Eq. (|4.8Q was derived. 
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The Fourier transform of the equal-time m — > ml particle Green's function: 

which we denote by G(E' ,p[,P2, ■ ■ ■ ,Pm'i EiPhP2, ■ ■ ■ ,Pm), is given by the formula: 

dz^dz^ ■ ■ ■ dzWdzW'dz®' ■ ■ ■ dzW 



G(E',p[,p2, . . . ,p^,; E,p u p 2 , . . . ,p m ) = J 



{2k) 



m+m' — 2 



m m 

*G{E' - z {i) ',p[, z {2) ',P2, • • • , z {m,) \p' m ,- z {i \pi, z (2 \p 2 , . . . , z {m \p m ), 

i=2 i=2 

where G(p' 1 ,p' 2 , ■ ■ ■ ,p' m >',Pi,P2, ■ ■ ■ ,Pm) is to be calculated using the Feynman rules in mo- 
mentum space for the m — > w! particle Green's function. These are the Feynman rules for 
the two-nucleon Green's function, with changes to Rules 1 and 2 to accommodate the differ- 
ent numbers of particles in the initial and final states. In Rule 1 we replace "two incoming 
and two outgoing" by "m incoming and m! outgoing" and Rule 2 becomes: 

"2. Assign to the incoming external particle legs four-momenta pi,p2, ■ ■ ■ ,Pm] assign 
to the outgoing legs four- momenta p'i,p 2 , ■ ■ ■ ,p' m >- Assign to the internal nucleon and pion 
legs momenta fci, . . . , ki, where I is the number of internal lines. Note that once the four- 
momenta of the particles in the initial state are chosen there are m'\ ways of choosing the 
four-momenta of the particles in the final state. All possible assignments must be included, 
as separate diagrams, with appropriate relative signs. If, for one such diagram, we make an 
assignment of the fermion four-momenta . . . ,p' n ), where rip is the number of fermions 
in the final state, then the relative sign of the other m'\ — 1 diagrams generated from it by 
permutation of the fermion four-momenta is (+1) for an even permutation of (p[, . . . ,p ; np ) 
and (-1) for an odd permutation." 

Therefore, the equal-time m — >• ml single-energy Green's function: 

G(E,p[,Pz, . . .,p^,,pi,p 2 , ...,p m ) 

may be found by calculating 



G(E' - zil) 'i P'n 2 {2) ', z {m,) ', p' m ,]E-Y, z {l) , Pi, z {2) z (m) , p m ) 

i=2 i=2 



according to these rules and integrating over z^ , z^ 3 \ . . . , z^ m \ z^\ . . . ,z^ m '^ . Again, this 
result is true regardless of whether or not the anti- nucleon approximation is used or not. 



VI. DERIVATION OF OLD-FASHIONED TIME-ORDERED PERTURBATION 
THEORY AND THE RESULT OF KVINIKHIDZE AND BLANKLEIDER FROM 

THIS PERTURBATION SCHEME 

In this section we will show how the above formalism allows the summation of all relative 
time-orders in an old-fashioned time-ordered perturbation theory diagram, and so derive the 
result of Kvinikhidze and Blankleider (KB) [|3T| . In the light of this calculation we will then 
comment on certain approximations made by KB in their work on the iViV — nNN problem 
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35fl . In order to achieve these goals it is necessary to first examine how to derive old- 
fashioned perturbation theory rules for the two-nucleon Green's function from the above 
Feynman rules. It will be clear from the structure of this derivation that we could generalize 
this argument to the m — > m! Green's function. 



A. Derivation of Feynman rules for time-ordered perturbation theory 

Consider the equal-time Green's function in the no-anti-nucleon approximation: 

G NN {i! ,x[,x' 2 \t,xx,x 2 ). (6.1) 

Define equal-time two-nucleon annihilation and creation operators in the Heisenberg repre- 
sentation via: 

ip N N{t,Xi,x 2 ) =4>{t,x 1 )$(t,x 2 ), (6.2) 
ip NN {t, Xi,x 2 ) = ip(t, x 2 )ip(t, xx). (6.3) 

The interaction and Schrodinger representation equal-time two-nucleon operators may then 
be found in the usual way. 

Currently, the perturbation expansion for G^N^t', x[, x' 2 \ t, x lt x 2 ) is written in terms of 
individual particle operators, i.e. if) T s, ijfs and 0fs. The first step in our derivation of 
time-ordered perturbation theory is to rewrite this perturbation expansion in terms of the 
two-nucleon operator, ippjpj, and the pionic operators, The only change which needs to 
be made in order to rewrite the perturbation expansion in this way is to prove that the 
nNN vertex: 

fi in t(x) = ig J d 4 x N d A x' N d 4 x 7r ifj(x' N )'y 5 Tip(x N ) ■ (£(xx)T(x - x' N ,x - x N ,x - x w ), (6.4) 

which we currently use, may be replaced by a vertex: 

?4rf(2)(z)= W J d 3 x'(J d A x N d 4 x' N d A x w tp NN (x° N \x' N ,x')>y 5 T 

x4>nn(x°ni %Ni %') ■ 0(x 7r )r(x — x' N , x — X N , X — Xn)) (6.5) 

in which not a single nucleon, but a nucleon pair, is destroyed and recreated. (See Figure 
for a diagrammatic representation of this new vertex.) The proof that, in the absence 
of anti-nucleons, rewriting the vertex in this way does not change the value of the Green's 
function is outlined in Appendix ||. This result then allows us to rewrite the expansion for 
Gnn^', %{,%2, t, Xi,x 2 ) obtained from Wick's theorem in terms of the free Green's functions: 

G i x N (t',x[,X2;t,x 1 ,x 2 ) = (0|T(^ 7V (t',f^,f2)^!v7v(^^i,^2))|0), (6.6) 
Gjji(y'iV) — (O|T(0y(2/')0i(2/))|O). (6.7) 

When we do this we obtain a set of Feynman Rules in coordinate space for this theory. 

Now we examine the consequences of making the no-anti-nucleon approximation. The 
first consequence of this approximation is that, similarly to the one-nucleon case, 
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G%(t', x[, x' 2 ; t, x h x 2 ) = if t > t'. (6.8) 

Since in the no-anti-nucleon approximation this Green's function must appear between any 
two vertices it follows that the time-order of all vertices is now determined. This is to be 
contrasted with the perturbation scheme used in the previous section where the relative 
time-order of most of the events occurring on different nucleons was not determined. In 
that perturbation scheme absorption and emission of exchanged pions were the only events 
occurring on different nucleons whose relative time-order was specified. In this perturbation 
scheme the relative time-order of all events is determined. 

Other than this change, however, the two perturbation schemes result in the same Green's 
function, since the rewriting of the equal-time two-nucleon Green's function in terms of free 
two-nucleon Green's functions does not change its value. Suppose then that we call the 
result of a particular diagram evaluated according to the rules of the previous Section: 

GNN(i){t j2?d t,xx, x 2 ), (6.9) 

where the subscript (1) indicates the use of one-nucleon Green's functions in the calculation. 
Suppose also that we call the result constructed by the rules obtained below for the theory 
with a two-nucleon vertex: 

GW ( 2)(f, x{, x 2 ; t, xx, x 2 ), (6.10) 

where the subscript (2) indicates the use of two-nucleon Green's functions in this calculation. 
Then it follows that: 

GnN(2){^> %'xi ^1) ^2) = GnN(1){^) x'xi ^1) ^2)5 (6-11) 

All relative TOs 

where "All relative TOs" indicates that the sum is over all possible relative time-orders of 
events on different nucleons, with the exception of the emission and absorption of exchanged 
pions, whose relative time-order is fixed as soon as the diagram to be considered is chosen. 
This result may also be deduced from Equation ( |B15|) which was derived in Appendix ||. 
Note that both Green's functions must be evaluated in the no-anti-nucleon approximation 
if this equation is to hold. We shall return to Eq. (|6.11| ) shortly. 

First, however, we must complete our development of time-ordered perturbation the- 
ory, for we have not yet explained how the rules usually used in the calculation of 
Gnn(2) {f, x[, x 2 , t, Xx, x 2 ) come about. 

Take the Fourier Transform of the Green's function: 

GNN(2)(t': ^2! Xx, X 2 ) , 

in order to obtain: 

G NN(2) (E' ',p[,p 2 ;E,px,p 2 ). 

Just as was done for Gjvjv(i) above, in Eq.( 4.14 ), energy-momentum conservation in this 
Green's function may be expressed via: 

G NN{2) {E',p[,p' 2 - E,px,p 2 ) = (2n) A 5(E' - E)5^ 3 \p[ + p' 2 - Pl - p 2 )G NN{2) (E,p[,p^p u p 2 ). 

(6.12) 

The rules for calculating GjqM(2){.E,p[,p 2 ,pi,p 2 ) are found from the rules in co-ordinate 
space by the standard procedure explained in Section |III[ The rules thus obtained are as 
follows: 
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Feynman Rules for the two-nucleon Green's function in terms of the free 
two-nucleon Green's function, in the no-anti-nucleon approximation. 

1. Draw all those topologically distinct, connected diagrams which have two incoming and 
two outgoing legs. Remember that different time-orders now contribute to different 
diagrams. 

In each diagram: 

2. Label the pion legs with momenta ki,...,kp, label the external nucleon pairs with 
combined energy-momenta (E,pi,p 2 ) and (£" ,p[,p 2 )- There are two possible ways of 
assigning the momenta p[ and p 2 in the final state, and both possibilities should be 
included as separate diagrams, with a relative minus sign. Label the internal nucleon 
pairs with combined energy-momenta: 

(E^, q f\ q f\ (EW, g f\ q f\ q^\ 

Assign isospin labels li, l 2 , ■ ■ ■ , l n to the n vertices. Note that because of our rewriting 
of 'Hint above, a new nucleon pair is regarded as being created whenever a vertex occurs 
on either nucleon. 

3. To the nucleon pair with combined energy-momentum (E, qi,q 2 ) assign a factor 

4. To the pion line with four-momentum k, joining vertices with isospin labels l e and l a , 
assign a factor GjfJ e (k). 

5. For each pion absorption vertex occurring on the nucleon with momentum qi insert a 
factor: 

(-tg)(2nY5(E' - k - E)5^(q> - k - - q^n) 

xT(E' - E N (q{), qi E - E N (%), q h k , k). 

For each emission vertex occurring under the same assumptions insert a similar factor, 
but with (k , k) replaced by (— k , —k). 

Here 

l= [l if< = 2. (6 - 13) 

Also, (E, qi,qi) is the combined energy-momentum of the nucleon-pair before the pion 
emission or absorption; (£", q(, qj) is the combined energy-momentum of the nucleon- 
pair afterwards, (k , k) is the four-momentum of the pion and / is the isospin label 
assigned to the vertex. 

6. Integrate over all internal energies and momenta: 

dEWd 3 q [ 1] d 3 q ? dE^d 3 q?d 3 q { 2 2) dE^-VdPqf-VcFqf-V d 4 h d 4 k 2 d A k r 



(2vr)7 (2tt)7 (2tt) 7 (2tt) 4 (2tt) 4 (2tt)^ 

and sum over all repeated isospin indices. 
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In these rules G^ N (E, q\, q 2 ) is defined by: 

(2tt) 7 5(E' - E)5 {3 \q[-q 1 )5 {3 \q' 2 -q 2 )G^ ) N {E,q l ,q 2 ) = J dtdt' d 3 Xl d 3 x 2 d 3 x[ d 3 x' 2 

xe^ E,t '-^<-^G%(t', x[, x' 2l t, x u x 2 )6(t' - t ) e -< m -^-^\ (6.14) 

When this Fourier transform is evaluated we find: 

77(0) - ^ im 2 Y, ai , a2 u ai {qi)u ai {qi)u az {q\)u a2 {q 2 ) ( . 

G ^ E ^ ~ E N {q l)EN {q 2 ) E+ - E^q,) - E N (q 2 ) " ( ° 5) 

At the same time G n f^(k) is defined exactly as in previous sections and so, in this no-anti- 



nucleon approximation, Gjf^{k) is given by Eq.( |3.17| ) 



Note that all the arguments in this section assume that the vertex V has no analytic 
structure in its energy variables. This assumption is justified since the presence of analytic 
structure would indicate that intermediate states in Y could be exposed. Because Y is the 
bare vertex for the theory it cannot have such intermediate states, and so it must have no 
analytic structure in its energy variables. 

This represents one way of formulating the Feynman Rules for the two-nucleon Green's 
function in this "time-ordered" or "old-fashioned" perturbation theory. 

There is another more usual way of writing these rules. Instead of representing a j'-pion 
intermediate state using one Green's function for the two nucleons and another j for the 
pions, we may write a single Green's function for all of the j + 2 particles in this state. 
This change can be accomplished either at the co-ordinate space level (in a similar but more 
complex way to the above change from an expression in terms of one-nucleon operators to 
one in terms of two-nucleon operators), or it can be done at the momentum-space level, as 
follows. 

Consider any diagram containing a j-pion intermediate state. Suppose that this diagram 
is evaluated using the above Feynman rules. We know from these rules that, if the isospin 
indices which are irrelevant to our argument are suppressed, the Green's function for the 
diagram may be written: 



G J (E',p[,p 2 ';E,p 1 ,p 2 ) = 

dEd 3 ptd 3 p 2 d 4 k w d*k^ ■ ■ ■ d A k^ 



F{E',p[,p' 2 - EjtJ,, fcW, fc^, . . . , k^)G%(E, 



PUP2) 



(2tt) 7 (2n) 4 i 

xG®(kM)G<?\k®) ■ ■ ■ Gf\kW)Fl{Ej x J 2 , fc« k^ 2 \ . . . , fctf>; E,p u p 2 ). (6.16) 

(See Figure |8|.) The Feynman Rules also ensure that energy is conserved in F and F'. 
In particular, we may extract a factor (2n)5(E — E — Y^ m =i &o"^) from F^ and a factor 
(2ir)8(E + Dm=i ^0 — E') from F. Furthermore, F and F^ have no poles in any of their j 
fco-variables, because the propagator for the mth pion, which is the only part of the expression 
to contain analytic structure in the variable k^, appears explicitly in our equation for G 3 
and so cannot form part of F and F*. 

Using this energy-conservation result and substituting for G^/ N and G^ gives: 
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G\E\p[,p'^E,p x ,p 2 ) 

, se/ „, ^ rd^p^k^d^-'-d^dk^dk^-'-dk^ 

= (2n)5(E - E) J — 

F(E', p{,p>; E-f: ki m) , 0i,02, , k^ , . . . , fcCi)) 

m=l 

im 2 E ai ,g 2 U ai (Pl)u ai (gO^gg (P2)^a 2 (P2) A 1 



E^E^) E+ - El=i K m) - Enfa) - E N (p 2 ) ^ 2^(A;M) fc M + _ Ww (£(m)) 

F'<(E-f2kt\0i,02,k^,k^,...,k^;E,p 1 ,p 2 ). (6.17) 

m=l 

Performing the &o integrations then gives: 

^ W ^ ^ „ _ _ ,„ j-,-, fd 3 p~ 1 d 3 p~ 2 d 3 kWd 3 kW---d 3 kW 
G J (E',p[,p^E, Pl ,p 2 ) = (2tt)5(E'-E) J (2^+6 

F(E>,p[,p>; E-f^ ^{k [m) ),0i,02, Ccv(£ (1) ), k^,co n (k^), k& , . . . , co n (k^), W) 

m=l 

im 2 1 Y, ai , a2 u ai (Pi)u ai (pi)u a2 (p 2 )u a2 (p 2 ) 



En{pi)En{P2) *=i 2w ff (fcM) £+ - E J m =i w*(tt m) ) ~ E N (fr) - E N (p 2 ) 
F\E - co 77 (k^),0 1 ,0 2 ; u;^ 2 )), fc< 2 >, . . . ,u«(k®), ^">; A, ft). (6.18) 

m=l 

Since this argument is valid for any j-pion intermediate state the above set of Feynman 
rules may be replaced by the following set of Feynman rules, which are those usually used 
for calculating the Green's function in old-fashioned time-ordered perturbation theory. 

Feynman rules for old-fashioned time-ordered perturbation theory 

1. As for Rule 1 above. 
In each diagram: 

2. Label the incoming (outgoing) nucleon pairs with momenta (pi,p 2 ) ((p[,P 2 )) and inter- 
nal nucleon pairs with momenta (q^, q 2 ^), (<fi 2 \ q^), ■ ■ ■ , (^l™ *\ 92™ including 
separate diagrams with a relative minus sign for the two possible momentum assign- 
ments in the final state. Label internal pion lines with momenta k^, k^, . . . , W \ 
Assign isospin labels h,l 2 , ■ ■ ■ ,l n to the vertices. Recall that a new nucleon pair is 
regarded as being created whenever a vertex occurs on either nucleon. 

3. If j pions are present at the same time as the nucleon pair, and they have momenta 

. . . , k^ while the nucleon pair has momenta q\ and q 2 , then associate, with all 
the particles present at that time, a Green's function: 



im 



2 



*T 11 - ; 7Z. 7T7^ — „ ,^ ^ "7*^7 IFwa, (6-19) 



ENiqjEnfa) 2u n (k( m y) E+ - E N {&) - E N {q 2 ) - El=i w v (kW) L\ 
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where l e and l a are the isospin indices of the vertex at which the Zth pion is emitted 
and absorbed. 

Here " j pions present at the same time" means that a vertical line drawn through the 
two nucleon lines intersects j pion lines. 

4. For each pion absorption vertex occurring on the nucleon with momentum q} at the 
same time as j spectator pions with momenta k±, . . . , kj are present, insert a factor: 



{-ig){2^f5^\q' t - k - ftX^Ws* - m )(2ic)* U 5^\k' m - kj^n 



3 

1 

m=l 



3 

m=l 



xT(E - J2 Mfy - E N ($),qt, E - E N (q\),u w (k),k). 



Here / is the isospin label assigned to the vertex. 

For an emission vertex occurring under the same conditions we insert a similar factor, 
but with k replaced by —k and 0J^{k) replaced by —u) n (k). 

5. Integrate over all internal momenta: 

fk^cPk^ fkp ggjV^ rf 3 g{ 2) rf 3 gf (Pqf-VcPqp-V 
(2tt) 3 (2^) 3 ' ' ' (2vr) 3 (2tt) 6 (2tt) 6 ' ' ' (2tt) 6 

and sum over all repeated isospin indices. 

B. The result of Kvinikhidze and Blankleider 

Having obtained the usual time-ordered perturbation theory rules for calculating 
Gnn(2)(E',p[,P2, E,pi,p 2 ) we may use Eq.( |6.11| ) to derive some of the results proven by 
KB using different means |3ll , |35| , |53| . 

Consider any time-ordered perturbation theory diagram, and suppose that the Green's 
function for the diagram is G/vjV(2)- Define a Green's function G NN ^ which is the sum 
of all Green's functions representing diagrams which differ from the diagram under con- 
sideration only in the relative time-order of processes occurring on different nucleons, but 
excluding those diagrams in which the order of emission and absorption of transmitted pions 
is different. I.e. define G NN ^ to be: 

G N N(2) = ^2 G NN ( 2 )- (6.20) 

All relative time-orders 



Once this definition is made we may take the Fourier transform of Eq.( |6.11 ) and use 
Eq.flOQ to obtain: 

G% N(2 )(E',P[,P2, E,P!,P2) = 



(27T) 
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where Gnn(i) must be evaluated in the no-anti-nucleon approximation. Or, if we remove 
the energy delta-function on both sides: 



G NN{2) (E,p[,p 2 ,p u p 2 ) = J -j—^G NN{1) (E - z',p[, z\p' 2 \ E - z,p!,z,p 2 ). (6.22) 

These are two of the crucial results of this paper. They state that the Green's function for a 
sum of time-ordered perturbation theory diagrams may be expressed as the integral over all 
initial and final relative energies of another Green's function, which is simply the covariant 
Green's function in the no-anti-nucleon approximation. This is the general result, now let 
us look at two specific examples. Firstly, we will consider an arbitrary disconnected diagram 
and show how Eq. ( |6.21| ) reduces to the original KB convolution integral result for this case. 
Secondly, we will examine one-pion exchange as an example of what happens when we apply 
Eq.( |6.22D to a connected diagram. 



1, Disconnected diagrams 

Consider an arbitrary disconnected diagram. By definition, no pion is transmitted from 
one nucleon to the other. Suppose that, in this diagram, processes PI occur on one nucleon 
and processes P2 occur on the other. Applying the Feynman Rules in Section |V| to this 
situation we see that: 

G^ N1 (E' - z',p[, z\p' 2 - E - z,p!, z,p 2 ) = (2tt) 8 S(z' - z)5(E' - E) 

x5^(p' 1 -p 1 )S^(p' 2 -p 2 )G P r 1 (E - (6.23) 

where G P T 1 (E—z,pi) and G P I 2 (z,p 2 ) are the single-nucleon Green's functions for the processes 
PI and P2 respectively, and are to be calculated from the rules in Section |TJ, applied in 
the no-anti-nucleon approximation. 

Compare this with the form of Gnn(2) f°r the same diagram: 

G^\ 2) (E',p[,p^ E,p h p 2 ) = (2n) 7 8(E' - E)5^(p' 1 - p 1 )8^{p' 2 - p 2 )G% N2 \E,n,p 2 ). 

(6.24) 

Applying Eq.( |6.22j ) to this specific case consequently yields: 

Glf m {EJ x ,p 2 ) = J ^G p N \E-z,p 1 )G p N \zJ 2 ), (6.25) 

which is KB's result up to a factor of i |JT|. This factor is missing because we are using 
propagators which differ from KB's propagators by a factor of i. If we replace all our 
propagators by KB propagators then, for an nth order diagram, a factor i n+1 appears on 
the left-hand side of Eq. ( |6.25[ ) and a factor i n+2 appears on the right-hand side. Therefore, 
we find that, for any diagram: 

Glt\ 2) {E,p u p 2 ) = ~J dzG p N \E - z^G^n), (6.26) 
where the propagators now agree with KB's propagators. 
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2. One-pion exchange 



Now we consider one-pion exchange. Suppose a pion is transmitted from one nucleon, 
which we label N2, to the other, which we label Nl. Suppose that before the pion is emitted 
from N2 (absorbed on Nl) certain processes PI (P2) take place on Nl (N2), and after the 
pion is emitted on N2 (absorbed on Nl) certain processes PI' (P2 r ) take place on Nl (N2). 
Then the Green's function constructed according to the rules in Section |^ is: 

Gnn(i)(E - z',p[, z' E - z,p u z,p 2 ) = (2n) 3 5 {i) (p 1 + p' 2 - p 1 -p 2 ) 
xG pl '(E - z',p[)(-ig)r(E - z',p[, E - z,p ljZ - z',p[ -px)G p \E - z,p x ) 
xG w (z - z',p 2 - p^G^ 2 ' (z' ,p 2 )(-ig)T(z' ' ,p' 2 ,z,p 2 ,z' - z,p 2 - p 2 )G^ 2 (z,p 2 ), (6.27) 

where all single-particle propagators G P T (po,p) are calculated using the rules giving the 
single-particle Green's function for a process P, given in Section |T|, applied in the no-anti- 
nucleon approximation, and we have suppressed all the spin and isospin structure of the 
result. 

Substituting this into Eq. (|6.22|) leads to: 

Gnn { 2) (#> PvP^PuA) = (2n) 3 5^ (p[ + p' 2 - Pl - p 2 ) J 

xG pl '(E - z',p[)(-ig)T(E - z',p[, E - z,p ± ,z - z',p[ - pi)G P1 (E - z,fr) 

xG w (z - z',p 2 -p>)G P2 \z',p 2 ')(-tg)r(z'.J>,z,p 2l z' - z,&- p 2 )G P2 (z,p 2 ). (6.28) 

If we again remove the relative factor of i between our propagators and those of KB we 
obtain, for an nth order diagram, a factor i n+2 on the right-hand side and i n on the left- 
hand side. Therefore, for any order diagram the result is: 

G E NN{2) {E,p[,p' 2 ,p l ,p 2 ) = 6 {3) {p' l +p' 2 -p 1 -p 2 ){-^-) 2 J dzdz' 

xGff'iE - z\p[)(-ig)T(E - z',p[,E- z,p ± ,z- z',p[ -pi)G^(P - z,p x ) 
xG^z - z',p 2 - p 2 )G^ 2 '(z' \p 2 )(-ig)T(z',po[, z,p 2 , z' - z : p 2 - p 2 )G^ 2 (z,p 2 ), (6.29) 

where, once again, the propagators now agree with KB's propagators. If one now sums this 
result over all possible processes PI, P2, PI' and P2' one obtains the sum of all time-orders 
of one-pion exchange from N2 to iVl with fully dressed propagators. 

G% N(2) (E,p[,p 2 ,pi,p 2 ) = 5 (3) (pi +p' 2 -pi -p 2 ){-^—) 2 I dzdz' 

xG N (E - z',p[)(-ig)T(E - z',p[,E- z,p 1: z- z',p[ -p x )G N {E- z,px) 
xG v (z- z',p 2 -p 2 )G N (z',p 2 )(-ig)r(z J ,p 2 ,z,p 2 ,z' - z,p' 2 - p 2 )G N (z,p 2 ), (6.30) 

where all of the propagators Gn are fully dressed. This is another result also obtained by 
KB, as an extension of their formula (|6.26|) [53]. As we shall see in Section VIII this formula 
for the one-pion exchange potential was also obtained by Klein and McCormick [4~3|. (See 
Figure |9| for a pictorial representation of this expression and the energy and momentum 
assignments involved.) 
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Clearly, this approach can be extended to two pion exchange, three pion exchange etc.. 
The advantage of such a method of calculating Green's functions is that it allows us to sum 
all relative time-orders in time-ordered perturbation theory by merely doing a convolution 
integral of the result obtained from the rules given in Section fV|, applied in the no-anti- 
nucleon approximation. 



C. The approximation of Kvinikhidze and Blankleider 

Equation ( |6.30| ) is the exact result for the sum of all relative time-orders for one-pion 
exchange from N2 to iVT, with fully dressed propagators. However, in their work on the 
NN — 7rNN problem KB use only an approximation to this result [35]. They represent the 
one-pion exchange Green's function by a product of two convolution integrals separated by 
an inverse nNN propagator: 

Gnn(2)( E 'Pi,P2iPuP2) ~ S {3) (p[ + p' 2 - p! - p 2 )[(-^-) 2 I dz'dz" 
y ' Ztxi J 

xG N (E - z\p[){-ig)Y{E - z' ,p[, E - z",p u z" - z',p 2 - p' 2 ) 
xG N (E - z\p l )G^z" - z',p[- p^G^z'^WG^N- 1 ^)}^-^-) 2 [ dz"'dz 

xGn{E — z,Pi)Gn(z, P2)], (6.31) 
where, once again, all propagators are fully dressed. KB have performed calculations which 



suggest that, at least for this case, the error in making such an approximation is small |53 
Be that as it may, we believe that such an approximation is physically flawed on at least two 
grounds, and consequently we expect it to produce invalid results, if not in the calculation 
of one-pion exchange, then in some other calculation. 



Firstly, to try and express one-pion exchange as in Eq. ( |6.3l[) violates conservation of 



energy in the following sense. If the energy arguments of the Green's functions are treated 
as off-shell energies of the particles involved, then the energy and momentum assignments in 
Eq. ( |6.31D correspond to the energies and momenta indicated in Figure ID[ These momenta 



are to be compared to those in the exact expression for the sum of time-ordered perturbation 
theory diagrams, depicted in Figure |S|. It is clear that in the exact expression the energy and 
momentum of each individual particle is conserved throughout the diagram, except at the 
vertices, where a pion interacts with a nucleon, thus modifying the energy and momentum 
of the nucleon with which it interacts. However, the tcNN vertices are constructed in such 
a way that energy and momentum are conserved there too; therefore, one can say that the 
energy and momentum of the individual particles is conserved in this exact expression for 
the sum of time-ordered perturbation theory diagrams. In the approximate expression, the 
momentum of each individual particle is still conserved, in the same sense. However, an 
inspection of the energy assignments reveals that, although the total energy is the same 
at all times during the pion exchange process, the energy of the three individual particles 
jumps suddenly and for no physical reason at the point where the two disconnected pieces 
are joined together to form a connected diagram. Since the relative energy of the particles 
is integrated over this jump in energy is possibly quite large, and so this violation of energy 
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conservation is bound to have an effect somewhere in any calculation performed using such 
an approximation. 

A second associated flaw in such an approximation lies in the diagrams it omits from 
the Green's function. The types of diagrams omitted are represented in Figure [0]. The 
justification for ignoring such diagrams is that, as mentioned above, KB have performed 



a calculation which shows their effect to be small [53]. This appears reasonable, until one 



carefully examines what physics is omitted and what is included in such an approximation. 
Observe that diagrams such as Figures [12] and [13] are still included in KB's calculation. 
Breaking down these two included diagrams and the omitted diagram into smaller pieces 
implies that KB's approximation amounts to always including the process on the left of 
Figure [14| but only sometimes including the process on the right. The process on the right 
of Figure O is excluded in some diagrams, such as Figure O, but included in others, such as 



Figure [13|. In other words, in certain diagrams KB ignore the Jennings mechanism |32|j33 



when a spectator pion is present. I.e., diagrams which contain the same physical processes 
are treated differently, with some time-orders included and some excluded. 

These two problems reflect a fundamental flaw in the KB approach. KB use a convolution 
integral in order to sum all relative time-orders in a disconnected diagram, an approach, 
which, as we have seen, leads to the calculation of a Green's function which is a no-anti- 
nucleon approximation to the full covariant propagator. In other words, in KB's model, 
disconnected pieces are calculated in an approach which is an approximation to covariant 
perturbation theory. However, when they come to construct convolution integral expressions 
for connected diagrams they do so by joining disconnected pieces, written as convolution 
integrals, together using an inverse nNN propagator, as in Equation ( |6.31| ). This method of 
joining disconnected diagrams together is consistent with time- ordered perturbation theory, 
not with covariant perturbation theory. Consequently, their technique for forming connected 
diagrams from disconnected pieces is inconsistent with their calculation of disconnected 
diagrams. This leads, firstly, to violations of energy conservation, since although the total 
energy E is conserved the relative energy is not constrained to be conserved, and secondly, 
to the unequal treatment of certain diagrams which only differ in the relative time-order 
of processes on different nucleons, which is, ironically, the problem the convolution integral 
was introduced to fix. 

If one wishes to join disconnected pieces together in a fashion consistent with the use 
of a convolution integral for disconnected diagrams one needs to reexamine the roots of 
the convolution integral technique, as we have done here. This leads to expressions such 
as Eq.( |6.29|) which are the exact sums of all relative time-orders in a set of time-ordered 
perturbation theory diagrams. But, such a result is only an equal-time no-anti-nucleon 
approximation to the full covariant perturbation theory expression. Therefore, we intend to 
persist with the covariant calculation in which all the physics is naturally included, rather 
than pursuing approximations with questionable physical motivation, such as those used by 
Kvinikhidze and Blankleider. 



VII. AMPLITUDES 

In this section we obtain the amplitudes which correspond to the two-nucleon Green's 
functions we have used so far. This is done via the standard procedure for obtaining ampli- 
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tudes from Green's functions: LSZ reduction [J36H , in which the external legs of the Green's 
function are amputated in order to obtain the corresponding amplitude. Again, it will be 
clear that the work of this section can be extended to any m — *> m! particle Green's function. 
For the two-nucleon system the S-matrix is defined to be: 

(pi,P2,out|pi,p 2 ,in) = (2tt) 8 5 (4) (p / 1 - pi)5 {A) (p 2 - p 2 ) 

+ (27T)^ 4 \p[+p' 2 -p 1 -p 2 )(p[,p' 2 \T NN \ Pl ,p 2 ), (7.1) 

where: 

(Pi,P2\Tnn\pi,P2) = G^ 1 (p\)G N 1 (p' 2 )G NN (p\,p 2 ]p 1 ,p2)G N l (pi)G N l (p 2 ), (7.2) 

and we have ignored spin and isospin indices. Note the use of dressed nucleon propaga- 
tors Gjv here.0 Firstly, we apply this formula to the general two-nucleon Green's function 
obtained from the rules in Section |TV] : 

G NN(1) (E - z',p[, z',p 2 ; E - z,p u z,p 2 ). (7.3) 

Note that we have defined a total energy E = pf + p 2 = p\ + p 2 and relative energies 
z' = p 2 ,z = p 2 in order to facilitate the derivation of a relation between the amplitude 
obtained from this Green's function and the one obtained from the time-ordered perturbation 
theory Green's function. When the free one-nucleon Green's function is substituted in, this 
equation becomes, in the no-anti-nucleon approximation: 

G NN{1) (E - z',p[, z',p 2 ; E - z,p u z,p 2 ) 



E+-z'-E N (p[)z' + -E N {p^) 

i i 
E+ - z - E^Sh) z+ - E N {p 2 



xT NN(l) (E - z',p[,z',p 2 ; E - z,p\, z,p 2 )— : j^ — ftt^T' ( 7 - 4 ) 



where E^ may contain terms representing the dressing of the nucleon. (The positive-energy 
projection operators may be omitted here, and throughout the rest of this section, since they 
have no effect, due to our use of the no-anti-nucleon approximation.) Secondly, consider any 
time-ordered perturbation theory Green's function: 

G NN{2) (E,p[,p 2 ,p u p 2 ) (7.5) 

LSZ reduction applied to G nn {2){E , p[, p 2 , pi, Pi) implies: 

G NN{2) (E,p[,p 2 ;pi,p 2 ) = 

1 

T NNm (E,p[,P2,Pl,p 2 )^- — T^TT ^ T^TT , (7.6) 



E+ - E N (p>) - E N {&) — ^'^^^^E+ - E N {&) - E N (p 2 ) ' 



3 Since we are using dressed nucleon propagators we should include factors of Z~z, where Z is 
the wave function renormalization, if we wish to get the LSZ reduction formula exactly correct. 
However, the inclusion of such factors makes no difference to the overall result and so we omit 
them in order to make the argument as clear as possible. 
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where, once again, En may contain terms representing the dressing of the nucleons.FI Above 
y 

we defined G NN ^ to be the sum of all Green's functions differing only in the relative time- 
order of processes which occur on different nucleons, excluding transmitted pion emission 
and absorption. Now we also define: 

T NN(2) = Tnn{2)- (7-7) 

All time-orders 

Since Eq.( |7.6|) holds between each term in this sum over all time-orders, it must relate 

Gnn(2) t° ^NN2 as wen - Therefore, Eq. ([7.6|) for G NN ^ and Eq.(|7.4|) for Gnn(i) ma Y be 

^^^^ y 

substituted into Eq.( J6.22| ) in order to obtain the relationship between T NN ^ and Ijvjv(i)- 

When this is done we find: 



~s „ _ _ f dz dz' 1 1 

T NN(2) (E, Pl ,p 2 , Pl , P2 ) = J 7^-^(^373^) + ^-E^ 

1 1 



xT NN{1) (E - z',p[, z',p- E - *,&>*,PM e+ _ z _ En(Pl) + j—Eim * (7 ' 8) 
If we wanted to calculate a physical quantity from a sum of time-ordered perturbation theory 

Y 

diagrams in which all relative time-orders are included, we would need to calculate T NN ^ 
on-shell. However, it may well be easier to calculate Ijvjv(i) and then use the relation between 

Y Y 

Tnn(i) an d T NN , 2 -> to obtain T NN , 2 ~. on-shell. 

In order to find out if this is feasible we set E = E on with 

E on = E N (p x ) + E N (p 2 ) = E N (p{) + E N {&). (7.9) 

For this value of E Eq.( [7.8|) reads: 

dzdz' ,. . 1 



TNN(2)(E m ,p[,P 2 ,Pl,P2) = J — |^disc( 



But: 



Therefore 



(_ 2 «)2 K z ' — E N (p 2 ) 

xT NNil) (E on -z',p[,z',pi;E on -z,p 1 ,z,p 2 )disc( ^-t^t). (7.10) 

z- E N (p 2 ) 



disc( ) = -2m5(z - w). (7.11) 

z — w 



T NN(2)( E ° n iPiiP2iPiiP2) = T NN(i)(E N (p[),p[,E N (p 2 ),p 2 ]E N (p 1 ),p 1 ,E N (p 2 ),p 2 ). (7.12) 

This is a more general version of a result obtained by Kvinikhidze and Blankleider by 
different means ||35|| . It states that, on-shell, the amplitude obtained from the no-anti-nucleon 
two-nucleon Green's function is equal to the sum over all time-orders of the amplitudes 
obtained from time-ordered perturbation theory. 



4 Note that the two-nucleon propagator used in this equation does not contain the full dressing 
on both nucleons, since certain time-orders are excluded from it. As discussed in Section |, this 
inadequate dressing was the original reason for KB's suggesting the use of a convolution integral 
for the free two-nucleon Green's function [31|. 
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VIII. THE METHOD OF KLEIN, LEVY, MACKE AND KADYSHEVSKY 

(KLMK) 



In the 1950s and 60s there was considerable interest in techniques for reducing the Bethe- 
Salpeter (BS) equation to a three-dimensional integral equation. In this section we explain 
how this work is connected to our work above and in so doing reveal that the KLMK method 
suffers from similar inconsistencies to those of the KB method discussed above. 

Initial work by Salpeter showed that if an instantaneous interaction was used in the 
Bethe-Salpeter equation a three-dimensional integral equation for the positive-energy com- 
ponent of the wave function was obtained ||37|| . This work was then extended by Levy and 



Klein, who showed that, in the ladder approximation, the BS equation could be approxi- 
mated by the following integral equation for the wave function (p ++ (p) fl38,40,13|: 



(E - 2E N (p)) (p++ (p) = A«(p)AV 2) (-p) / d 3 kV OPE (E,p,p - k)^ ++ (p- k), (8.1) 
where the wave function ip ++ {p) is defined by: 
(E-2E N (p))<p ++ (p) = 

(-2 m )^E + Po - E N (p))( l -E-p - E N {-p))A^(p)Af(-p)i;(p), (8.2) 

with p the relative four-momentum and E the total energy in the center-of-mass frame. The 
kernel Vope(E,p,p — k) is given by: 

1 ,o r „ . 1 



Vop E {E,p,p-k) = -\{E -2E N {p))[{—) 2 J d 3 kdk dp n 



E + p -E N (p) 



2' 

X -j— L^_T (1) (p p _ k ) —^T^ (-p, - P + k) 

±E - p - E N (-p) k 2 + /j, 2 

x- ^- =, }{E-2E N (p-k)). (8.3) 

t\E + p - k - E N (p- k) jE - p + k - E N (k - p) 

If we regard E^ as containing a self-energy term for the nucleon then this kernel contains 
both fully dressed vertices and fully dressed propagators and we may write: 

V OPE (E,p,p - k) = {E- 2E N (p))Gl N{2) (E,p, -p,p- k,k- p){E - 2E N (p- k)), (8.4) 

where G^ N , 2 JE,p, —p,p — k,k — p) is the sum of the convolution integral Green's functions 
for pion exchange from Nl to N2 and N2 to Nl, which are given respectively, by Eq.( p.30| ) 



as written, and by the same equation with Nl and N2 swapped. Note that in order to 
establish this fact we need to perform the following substitutions, which merely correspond 
to changing the frame of reference to the centre-of-mass frame for the two-nucleon system: 

p = p[ = -p 2 "; (8.5) 

Po = \e- z'- (8.6) 

k = p2-p' % =p[-p\\ (8.7) 
k = z - z'. (8.8) 
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Using Eq.( [7.6|) we see that this Vope is the amplitude for the sum of all the time-ordered 
perturbation theory diagrams which represent the different time-orders of fully-dressed one- 
pion exchange. 

Klein went on to generalize this result, showing how the BS equation: 

M — I + IGM, (8.9) 

where M and / are regarded as operators in the Hilbert space of energy-momentum states 
\p) could always be written as a three-dimensional equation of the form: 

T{E) = V{E) + V(E)G++(E)T(E), (8.10) 

where 

G+ + N (E) = A i + ) ^+ ) G NN (E), (8.11) 
with the t-matrix T(E) defined by: 

T(E) = (^) 2 G N UE)[[dp' d Po G^\^E-p' )G^\]-E + p' ) 

xM(p' , Po ; E)G { ^(^E - Po)G { +\±E + p )]G^ N (E), (8.12) 

and the potential V(E) given by: 

V(E) = V\(E) + V 2 (E) - V 1 (E)G++(E)V 1 (E) + . . . , (8.13) 

where: 

V,{E) = (^-) 2g nn(E)[J dp' Q dp,G^\\E - p' Q )G { +\\E + p' ) 

xI(p' ,p 0] E)G^\ 1 -E -p,)G { +\ l -E + Po )}G N 1 N (E), (8.14) 

V 2 (E) = (JLy»G]&(E)[/ dp' dk dp G^ (^E - p' Q )G { +\\E + p' Q )I{p'„k -E) 

xG(k , E)I(k ,p ; E)G { +\^E - p )G^\^E + p )]G N 1 N (E) ■■■■ (8.15) 

All operators in these formulae are now considered to be operators in three-dimensional 
momentum space. 

Note that the propagator Gnn(E) is the time-ordered perturbation theory propagator: 

Gnn(E) = — — | — . (8.16) 

If dressed particles are involved the use of this time-ordered perturbation theory propagator 
is open to question a la KB | 3T| . But, we ignore this difficulty for the present and assume 
that this Gnn(E) is the correct propagator to use, since any such problem may be fixed by 
making a minor modification to the approach outlined here. 

Now if the interaction kernel / contains only the single pion exchange that gives the 
ladder BS equation then: 
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{p'\yi{E)\p) = V OPE {E,p\p). (8.17) 

This result was also obtained by Kadyshevsky and his collaborators [ 47 J48| , by proceeding 
from the Lippman-Schwinger equation for the amplitude. Klein's early work focused on 
deriving this result for the wave function BS equation, although he later also derived it for 
the BS equation for amplitudes [Q. 

Klein also explained how to obtain a time-ordered perturbation theory expression for 
V(E) from the results given here PJ. This explanation indicated that the potential V(E) 
could be regarded as the sum of all possible time-ordered perturbation theory diagrams for 
the process under consideration. That is to say, given any BS equation interaction kernel 
/, there are two equivalent ways of obtaining the three-dimensional integral equation (|8.10|) 
corresponding to the BS equation (|8.9|): 



1. Use Eq.( fg.l3j ) in order to derive the equivalent energy-dependent potential V(E); 

2. Sum all two-particle irreducible time-ordered perturbation theory diagrams allowed 
for this /; the result is V(E). 

That these two methods are both equivalent to the BS equation is a manifestation of the 
result derived above in Sections [TVHVTl: both the convolution integral and time-ordered 
perturbation theory may be obtained from field theory in the equal-time no-anti-nucleon 
approximation. 

Although in theory both of these two methods provide three-dimensional integral equa- 
tions containing exactly the same information as the BS equation, in practice it is not 
feasible to obtain the full V(E) which is necessary in order to make the equivalence exact. 
Consequently approximations to the potential, such as that defined by Equation ( |3.3|) , must 
be made. Sometimes additional modifications to the propagator are also made in order to 



further facilitate the calculation, see e.g. [48]. However, here we merely wish to examine the 
quality of the approximation: 

V(E) = V 1 (E), (8.18) 

which is obtained by taking the first term in the series (|S.13|) for V when the BS equation 
is being used in the ladder approximation, and so I consists only of one-pion exchange. If 
dressed propagators are being used the time-ordered perturbation theory diagrams which 
are omitted from this approximation but included in the ladder BS equation include those 



shown in Figure 15. These omitted diagrams are to be contrasted with those in Figure 16 



which are included in the calculation using Eqs. flS.lOD and ( |8.18[ ). It is clear that, as is the 



case in the KB method and the current NN — nNN equations, different time-orders of the 
same physical process are treated differently. 

Connected to this problem is the problem of relative energy non-conservation, which 
occurs here in a similar way to that in which it occurred in the KB approximation discussed 
in Section |V1 C| . If we examine the energy assignments for the two-pion exchange diagram 
generated by Equations ( 8.10|) and ( 8.18|) we find that they correspond to those drawn in 



Figure [T7|. Although the energy assignments for one-pion exchange are now correct (compare 



with Figures |9] and [10]) there is still a sudden jump in energy where one one-pion exchange 



diagram, calculated using a certain pair of relative energies, is joined to another one-pion 
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exchange diagram, which is calculated using different relative energies. Again, as discussed 
above for the KB approximation, this jump, could, in principle, be infinite, since the relative 
energy variables k Q and k' are integrated over. Consequently, an examination of KLMK's 
techniques shows that although their work contains many of the same elements as those 
pursued in this paper, ultimately it has a different purpose, since it places a higher value 
on obtaining a tractable three-dimensional integral equation than on including all possible 
time-orders of the relevant processes. In this sense it is more akin to the work of KB than to 
the approach espoused here, since both KB and KLMK restrict themselves to considering 
only three-dimensional integral equations and using energy-integrations to sum certain pieces 
in their essentially three-dimensional theories. They pay for this restriction however, since 
both approaches violate relative energy conservation and neither includes all time-orders 
of physically relevant processes. If these two problems are to be eliminated a fully four- 
dimensional calculation must be pursued. 



IX. CONCLUSION 

Current models of the NN — ttNN system are unsatisfactory in two ways. Firstly, 
they fail to correctly predict the experimental data. Secondly, because all calculations are 
done, in some sense or another, in a time-ordered framework, diagrams which are merely 
different time-orders of the same process are treated unequally — the most famous example of 
this being the non-inclusion of the Jennings diagram in the current NN — wNN equations 
|32] , |33|| . It is to be hoped that a resolution of the second problem will, by producing a 
theoretically consistent model, lead to a diminution of the first problem. Therefore, the 
question is how to consistently include all time-orders in a calculation. 

In this paper we displayed a scheme for doing just this. Starting from a Lagrangian 
closely related to those of Quantum Hadrodynamics |7],[8| we explained how to calculate 
fully covariant one-nucleon, one-pion, two-nucleon and m — > m'-particle Green's functions, 
either with or without anti-nucleonic degrees of freedom. It was then shown that, in the 
absence of anti-nucleonic degrees of freedom, any covariant perturbation theory diagram 
may be expressed as the sum of a set of time-ordered perturbation theory diagrams, pro- 
vided we restrict all particles to have the same time in the initial and final states. This 
result made it clear that Kvinikhidze and Blankleider's original convolution integral formula 



for disconnected diagrams |31| is merely a manifestation of the more general relationship 
between covariant perturbation theory and time-ordered perturbation theory. This fact al- 
lowed us to point out the inconsistency in the KB approach to the NN — irNN problem 



3"5] , which leads to a violation of relative energy conservation. To complete our derivation of 



the properties of convolution integral theory we explained how LSZ reduction |36| could be 
performed on convolution integral Green's functions. We then observed that KLMK's work 
[p8|--[48| , which also used a convolution integral to calulcate the one-pion exchange potential, 
has problems similar to those which occur in the KB approximation. These problems occur 
because any approximate three-dimensional integral equation, which involves an energy in- 
tegration for the propagators, as in the work of KB, or for one-pion exchange, as in the work 
of KLMK, results in relative energy non-conservation. This is connected to the fact that 
such three-dimensional integral equations always omit certain time-orders of some physical 
process for which other time-orders are included. And, if there is one thing the failure of 
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the current NN — nNN equations teaches us, it is that we treat different time-orders of the 
same physical process differently at our peril. As stressed above, the only way to circumvent 
this difficulty and consistently include all time-orders of the relevant physical processes in 
a natural way is to use time-dependent perturbation theory and derive four- dimensional 
integral equations. 

Having developed these points in this work the next task is the application of the per- 
turbation theory developed in the first five sections of this paper in order to derive fully 
covariant scattering equations for systems of nucleons and pions. The technique chosen 
to do this is the classification-of-diagrams method of Taylor |15[| . Because there are some 
subtleties in the application of Taylor's method to perturbation schemes other than time- 
ordered perturbation theory we will explain Taylor's technique fully and elucidate these 
subtleties in a forthcoming paper [ 53J . In a third paper the results developed in this second 



paper will be applied to the iViV — ttNN system [55 
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APPENDIX A: CALCULATION OF THE ONE-PARTICLE FREE GREEN'S 

FUNCTIONS 

In this appendix we calculate the one-nucleon and one-pion free Green's functions. We 
then apply the no-anti-nucleon approximation to these Green's functions, and hence obtain 
the form of the free Green's functions in the no-anti-nucleon approximation. 



1. Calculation of the full covariant Green's functions 

We begin by noting that the Green's functions in question are defined by: 

G i °\x',x) = (O\T(ij I (x'0(x))\O), (Al) 
Gjfi (x',x) = (O|T(0 J / (x / )0[(x))|O). (A2) 

Firstly, consider the one-nucleon Green's function. Using the definition of the time-ordering 
operator T, the expansions of the Schrodinger representation field operators i/j(x),ip(x): 

^ = I yf^^TA EM^KWe** + dl(k)v a (k)e- rk % (A3) 
J (27TT E N (k) a 

= / T^^TA Y.iblm a (k)e- rk -* + d a {k)v a {ky % -% (A4) 

J (27Tf E N {k) a 
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and the connection between the Schrodinger and Interaction representations gives: 

G#V> x) = G% )+ {x\ x)9(x' - x ) + G { ° r (x', x)9(x - x' ); (A5) 

with: 



J (An) hj N {p)h N {p ) 

( )- { , f d S pd 3 p' m 2 Eqq' Va(p)Vg'(p' 



(2tt) 6 E N {p)E N {p>) 
<n *)=/ ^ ^ f;g ) f^l P ' ) (0\d a {p)e^^dim0)e i ^^, (AT) 

where we have simplified G$ and by eliminating those terms which are zero due to: 



6«(p)|0> = 0, (A8) 
d a (p)|0) = 0. (A9) 

Now consider G^\p), defined by: 

CSJV.P) = (27r) 4 ^ 4 )(p-p')G'g ) (p) (A10) 

where: 

Gn\p',p) = J d i xd i x'e^'G^\x',x)e-^ = G% )+ {p' ,p) + Gf>',P)- (All) 

Here G$ + (p',p) and are the Fourier transforms of G$ + (x', x)Q(x' Q — x ) and 

G$ (x',x)6(x — x' Q ) respectively. Evaluation of G$ (p',p) gives: 



G { n + (p',p) = J dx dx' 



,™ Ea,a' U a '(p')u a (p) 



E N {p)E N {p>) 

(0\b a/ (p')e- iHK( ^- Xo) bl(p)\0)9(x' o - xo)e i(p '° x '°- poXo \ (A12) 



Transforming to relative and average times: 



T = ^±^,r = x' -x 0] (A13) 



and noting that dx dx' = dTdr gives: 



= / a**-** f" ir ^ (0MpV«^>-"*M(p)|0> (A14) 

J jo En{P )En{P ) 

= ^^-V ^ 2 ^^^^ M0\aApy^- HK)T aUp)\0) (A15) 
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where Hk\p'ch') = E^{p')\p'a') . Now the Fock space states \pa) are normalized such that: 

E N (p) 



(p'a'\pa) = 5 aa , (2irf ^±5^ (p'-p), 



ni 



hence: 



gT{p\p) = {2-kY^\p'-p) 



im Y,a.Ug{p)u a {p) 

E N (p) p+-E N (p) 



(A17) 



(A18) 



Therefore, 



GP ' (p) 



im T. a u a (p)u a (p) 
E N (p) p+-E N {p) ' 



(A19) 



in agreement with Garcilazo up to a factor of i ||56|| . This result is also in agreement with 
van Faassen, again up to a factor of i. van Faassen absorbs the factor of into his spinor 
normalization [57]. 

A similar calculation for G$ (p) yields: 



im T, a Va(-p)Va(-p) 

E N (-p) Po+E N (-p) 



(A20) 



again in agreement with Garcilazo. Noting that En(—p) = E^{p) these two results may be 
combined in order to produce: 



im 



E N (p) 



Eg U a (p)u a {p) Eg Vg(-P)va(-P. 

Po ~ E N (p) 



Po +E N (p) 



(A21) 



which may then be simplified to yield: 



p 1 — m 



(A22) 



in accordance with Eq.( |3.12|) . 

We now turn to the free pion Green's function: 

G^ ) (x' ) x) = (O|T(0j(x')^(x))|O). 

Using the expansion of the Schrodinger operator: 

d 3 k 



(x 



ai (k)e ik -* + a\(k)e- lk - x ) 



(A23) 



(A24) 



and the relation of the interaction and Schrodinger operators gives, by a similar argument 
to that for the free one-nucleon Green's function: 



Gj$(xf, x) — Gju (x',x)9(x' — x ) + G w fi (x',x)9(x - x' ), 



'TTJl 



(A25) 



with: 
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G$ + W,x) = J ^ -J^ p (0|a J -(^) e -^K-»o) a t ( g ) | 0)e ^-^-^ > (A26) 

G$~V,x) = I J- (0|a,(fc)e-^^-o)a](fcO|0)e- ( ^ fc "" ,) . (A27) 

If we define (A;', fc) and (A;', fc) to be the Fourier transforms of Gj^j (x',x)9(x' — 

Xq) and (x',x)9(x — x' ), we may obtain, by a similar procedure to that used for the 

free one-nucleon Green's function above: 

G^\k', k) = (2n)W\k' - *)—*-=- ^ (A28) 
Gjji~ {k 1 , k) = - {2*)W(k> - k)-^- - (A29) 



Note that once again, these formulae agree with Garcilazo up to a factor of % |56j . Note 
also that: 

G$~{k\ k) = G^i-k, -k'). (A30) 

If we use the fact u^ik) = uv(— k), these may be combined to give the total Green's 
function: 

G${V, k) = ^vOV 4 ^' - k)Gjf>{k). (A31) 

with: 

G -*M = p^4- < A32 > 

in agreement with Eq. (|3.13| ). 



2. The no-anti-nucleon approximation 



So far we have not made any approximations in the calculation of these Green's functions. 
In this section we examine what happens if we make the no-anti-nucleon approximation, i.e. 
take the limit v(p) —>■ and v(p) — ► 0. Examining Eq. (|A7|) reveals, in this limit, we obtain: 

G < § r {x',x)=0, (A33) 

and so, 

G$(x',x) = G < S )+ (x',x)9(x' - x Q ), (A34) 

in the no-anti-nucleon approximation. This proves Eq. (|3.15|) . 
Furthermore, since: 

GgV.P) = G^\p\p) + gP~(p', P ), (A35) 
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and G$ (p,p') is zero in this approximation, it follows that: 

G { »\p\p) = G% )+ {p\p). (A36) 
Consequently, Eg. ( |A19| ) implies that, in the no-anti-nucleon approximation: 



G%> ( p ) = ^r ^T^ , ( A37) 

NKFJ E N (p) pt~E N (p) 1 ; 



as claimed in Eq. fl3.16l) . 



APPENDIX B: REWRITING THE PERTURBATION EXPANSION OF THE 
TWO-NUCLEON GREEN'S FUNCTION IN ORDER TO OBTAIN 
OLD-FASHIONED PERTURBATION THEORY 

In this Appendix we establish that the vertex containing single-particle operators, i.e.: 

'Hint(x) =ig J d A x N d 4 x' N d 4 x n t/j(x' N )'y 5 Til'(x N ) ■ c£(x n )T(x - x' N ,x - x N ,x - x n ), (Bl) 

may be re-expressed in terms of two-nucleon operators, which in the Schrodinger represen- 
tation are defined by: 

iI>nn(x 1 ,x 2 ) = iI>(x 1 )iI)(x2), (B2) 
ip NN (xi,x 2 ) = ^(^M^i)- (B3) 

It immediately follows that the equal-time two-nucleon Green's function may be rewritten in 
terms of these two-particle operators without changing its value, as was claimed in Section 

o 

Claim: Suppose that the state {S) 1 contains two nucleons, and both these nucleons in 
IS) 1 are created at the same time. Define 'H\ nt {2) by: 

nl m (y , y, y') = ^\yo, y')nUy)^(y , y') (B4) 
= ig J d^yNd^y'^y^xiy'x^y^y'hbT 

x^N N (yNo,VN:y') ■ ¥ r d/ 7 r)r(y - y' N ,y - y N ,y - y w ). (B5) 

Then: 

nUvW = J d'y'Hint^M'W. (B6) 
(Naturally these statements, if true in one representation, are true in any representation. 



We have chosen the interaction representation, because, as will be seen in Eg. ( \B14) , it is 
operators in that representation which the Green's function is expressed in terms of.) 

Proof: 

Consider: 
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/ dZy'Hintwiyo^y'W, (B7) 

where \S) 1 is a state in the interaction representation containing an, as yet, unspecified 
number of nucleons, all of which were created at the same time. Firstly we note that, using 
the anti-commutation relations of the ips and ^s, and the definition of TC( nt , we may obtain: 

n inm (y , y, y') = HUvo, y)e iH «^{y')^{y')e- iH ™ - HUvo, y)S {3) (y - y'). (B8) 

Now if the state \S) 1 consists of nucleons which are all created at the same time t and pions 
created at unspecified times it may be written, in the interaction representation, as: 



l^) 7 = II ^ f7 (^) II e^-V(^)|0), (B9) 

h=l j=l 

where and n s are the numbers of nucleons and pions in the state \S) and 
(t, Xi), . . . , (t, x u n) and z\, . . . , z n \ are the space-time points at which these nucleons and 

s s 

pions are created. Clearly, the operator (/ d 3 y'ip(y')ij(y')) commutes with all the pionic 
operators, and, it may also be shown to commute with Hk, leaving us with: 

j3 / 



d y'Hint(2)(yo,y,y')\S) = 

nUvo,y) II <f\^V Hk \ I W)^(f')) n^)l°> -nLM)\s)- ( B1 °) 

k=l J j=l 

Furthermore, from the anti-commutation relations of the ips and -0s we find that: 

' d 3 yW)^(y')Mx) = ^{x)+^{x){J d 3 y^{y')^{y')). (Bll) 
This implies that: 



AT N N 

l s s a s 



d 3 y^(y')^(y')) ]J ^(xj) = n» J[ ^(xj) + J d 3 y^{y')^{y')), (B12) 

3=1 3=1 3=1 

suggesting that (/ d 3 y'ip(y')if)(y')) may be interpreted as the fermionic number operator, as 
we would expect. Substituting this in Eq. (|B10| ) and simplifying the result gives: 

J d 3 y'lii ntm {y ,y,y')\Sy = {n N s - WUvoMS) 1 , (B13) 

which, if rig = 2 proves Eq. (|BT)|) , thus proving the claim. 

Now consider the definition of Gm^{x'^ x 'i'i x i-> x 2)-> Eq.( |4.1| ). By methods covered in most 
text books on field theory, the following expression for Gnn(xi,x 2 ; x[, x' 2 ) may be derived: 

00 (—i) n r 

G NN (x[, x' 2 \ x 1 ,x 2 ) = V j— / d^y x d^y 2 • • ■ d 4 y n 

x(0|T(^(x' 1 )V / (4)^L(^) • • .HL(2/2)HL(2/i)? 7 (^i)^(^))|0) co „, (B14) 
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where the con indicates that those contributions to the matrix element which contain 
vacuum-vacuum subdiagrams are to be ignored. Suppose that we make the no-anti-nucleon 
approximation. Suppose that we also consider only the equal-time Green's function, i.e. we 
set x\ = x\ = t and x\ = x® = t'. We begin by using the definition of T in order to 
rewrite this Green's function as a sum over all possible time-orders of the field operators in 
the vacuum expectation value. Each time-order is written as a separate term in the sum, 
with the integration over y®, y®, . . . , y® being suitably restricted. 

We wish to prove that, if we consider any term in the sum over all time-orders, replac- 
ing Hi n t(y) by / d 3 y r Tli n t(2)(yo,y,y r ) does not change the value of that particular term. So 
consider any one term in the sum over all time-orders. Firstly note that if the presence of 
anti-nucleons is forbidden then the time t must be the earliest time, and consequently the 
operators ip {x\) and ip (x%) must be the first to act. Now the state if) (x\)ip (x2)|0) has a 
nucleon number of two, and since nucleon number is conserved in this theory, it follows that 
the number of nucleons in the state on which Timt is acting is always two. Furthermore, if all 
vertices 7i% n t{y) U P to the jth one have been replaced by vertices / d 3 y'Hi n t(2)(yo, y, y') then 
it is clear that the jth interaction Hamiltonian Hi n t{yj) acts on a state in which two nucleons, 
created at the same time, are present. Consequently, TCi n t(yj) may also be replaced by a ver- 
tex J d^y'pi.int^iy^yj-iy'j)- Since the first TCi nt (y) acts directly on the state tp (xi)tp (^2)|0) 
which clearly contains two nucleons, which, by assumption, are created at the same time, 
it follows, by induction, that replacing all vertices Hi n t{y) by vertices J d 3 y'Hi n t(2)(yo,y,y') 
does not change the value of this term in the Green's function. Consequently, Eq.( |B14|) may 
be rewritten: 

G NN {x[,x' 2 ; xx, x 2 ) = — p E / dyld*y l d i y' l dyld z y2d?y'2 . . . dy° n d 3 y n d 3 y' n 

n=o n - AUTOs n 

X WNN{t,x[,X2)ni nt{2) {y°,y n ,yn) ■ ■ - ^L^^y^yi^y^NNiXl, X!, X2)\0) con , (B15) 

where Q is a region of integration appropriate for the particular time-order under considera- 
tion. Thus, we have proved that, in the absence of anti-nucleons, the equal-time two-nucleon 
Green's function may be rewritten in terms of the two-particle Green's functions without 
changing its value. 
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FIGURES 



FIG. 1. The so-called Jennings term, which is not included in the NN — ttNN equations. 



FIG. 2. The diagram on the left is not included in the NN propagator used in the NN — ttNN 
equations, while the diagram on the right is. 



FIG. 3. This diagram is currently included in the NN — ttNN equations, even though it is 
merely a different time-order of the excluded Jennings term. 



FIG. 4. The diagram on the right is not included in the Stingl and Stelbovics model, since it 
involves a three-pion intermediate state. This occurs even though it is merely a different time order 
of the left-hand diagram, which is included in the model. 



42 



FIG. 5. The ttNN vertex we will use in this work. 



FIG. 6. The diagram on the left contains a vacuum- vacuum subdiagram and consequently is not 
included in the perturbation expansion for the two-nucleon Green's function. However, the diagram 
on the right, although it is disconnected like the left-hand diagram, contains no vacuum- vacuum 
subdiagram and consequently is included in the perturbation expansion. 



FIG. 7. A diagrammatic representation of the two-nucleon vertex, in which one nucleon is a 
spectator. To be compared to the vertex in Figure || which is only a one-nucleon vertex. 



FIG. 8. A diagrammatic representation of the Green's function for a process with a j-pion 
intermediate state. 
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FIG. 9. The energy and momentum assignments in the exact convolution integral result for 
one-pion exchange. 



FIG. 10. The energy and momentum assignments in the Kvinikhidze and Blankleider (KB) 
approximation to the exact convolution integral expression for one-pion exchange, which is shown 
in Figure |9[ 



FIG. 11. An example of the diagrams left out of the calculation in KB's approximation to 
one-pion exchange. 



FIG. 12. An example of the diagrams included in the calculation in KB's approximation to 
one-pion exchange. 
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FIG. 13. A diagram containing the Jennings mechanism with spectator pion which is included 
in KB's calculation of one-pion exchange. 



FIG. 14. The process on the left is always included in KB's approximation to the exact convo- 
lution integral one-pion exchange formula, while the process on the right is not always included; it 
is forbidden in certain diagrams, e.g. Figure 11, and allowed in others, e.g. Figure 13. 



FIG. 15. These diagrams are an example of those omitted from the calculation of two-pion 
exchange in the KLMK method. Note that they are merely a different time-order of other included 
diagrams, see e.g. Figure 16. 
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FIG. 16. These diagrams are included in the calculation of two-pion exchange in the KLMK 
method. 



FIG. 17. The energy assignments in the KLMK method calculation of two-pion exchange. 
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